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Abstract 

In the simplest Higgs phase of gravity called ghost condensation, an ac- 
celerating universe with a phantom era (w < —1) can be realized without 
ghost or any other instabilities. In this paper we show how to reconstruct 
the potential in the Higgs sector Lagrangian from a given cosmological his- 
tory (H(t), p(t)). This in principle allows us to constrain the potential by 
geometrical information of the universe such as supernova distance-redshift 
relation. We also derive the evolution equation for cosmological perturba- 
tions in the Higgs phase of gravity by employing a systematic low energy 
expansion. This formalism is expected to be useful to test the theory by dy- 
namical information of large scale structure in the universe such as cosmic 
microwave background anisotropy, weak gravitational lensing and galaxy 
clustering. 



1 Introduction 



Acceleration of the cosmic expansion today is one of the greatest mysteries in both 
cosmology and fundamental physics Assuming that Einstein's general relativity 

is the genuine description of gravity all the way up to cosmological distance and time 
scales, the so called concordance cosmological model requires that about 70% of our 
universe should be some sort of energy with negative pressure, called dark energy. 
However, since the nature of gravity at cosmological scales has never been probed 
directly, we do not know whether the general relativity is really correct at such infrared 
(IR) scales. Therefore, it seems natural to consider modification of general relativity 
in IR as an alternative to dark energy. Dark energy, IR modification of gravity 
and their combination should be tested and distinguished by future observations and 
experiments [3 El 13 13 E| • 

From the theoretical point of view, however, IR modification of general relativ- 
ity H U EDI HD IH 113 13 113 CHI EE3 is not an easy subject. Most of the previous 
proposals are one way or another scalar-tensor theories of gravity, and are strongly 
constrained by e.g. solar system experiments jTSJ and the theoretical requirement 
that ghosts be absent ^3 123 121] • The massive gravity theory [22j and the Dvali- 
Gabadadze-Porrati (DGP) brane model [23 are much more interesting IR modifica- 
tion of gravity, but they are known to have macroscopic UV scales [23123- A UV scale 
of a theory is the scale at which the theory breaks down and loses its predictability. 
For example, the UV scale of the 4D general relativity is the Planck scale, at which 
quantum gravity effects are believed to become important. Since the Planck scale is 
microscopic, the general relativity maintains its predictability at essentially all scales 
we can directly probe. On the other hand, in the massive gravity theory and the 
DGP brane model, the UV scale is macroscopic. For example, if the scale of IR mod- 
ification is the Hubble scale today or longer then the UV scale would be ~ 1, 000km 
or longer. At the UV scale an extra degree of freedom, which is coupled to matter, 
becomes strongly coupled and its quantum effects cannot be ignored. This itself does 
not immediately exclude those theories, but means that we need UV completion in 
order to predict what we think we know about gravity within ~ 1, 000km. Since this 
issue is originated from the IR modification and the extra degree of freedom cannot 
be decoupled from matter, it is not clear whether the physics in IR is insensitive to 
unknown properties of the UV completion. In particular, there is no guarantee that 
properties of the IR modification of gravity will persist even qualitatively when the 
theories are UV completed in a way that they give correct predictions about gravity 
at scales between ~ 1, 000km and ~ 0.1mm. 
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Ghost condensation is an analogue of the Higgs mechanism in general relativity 
and modifies gravity in IR in a way that avoids the macroscopic UV scale [26J 1 . In 
ghost condensation the theory is expanded around a background without ghost and 
the low energy effective theory has a universal structure determined solely by the 
symmetry breaking pattern. While the Higgs mechanism in a gauge theory sponta- 
neously breaks gauge symmetry, the ghost condensation spontaneously breaks a part 
of Lorentz symmetry 2 since this is the symmetry relevant to gravity. In a gauge the- 
ory the Higgs mechanism makes it possible to give a mass term to the gauge boson 
and to modify the force law in a theoretically controllable way. Similarly, the ghost 
condensation gives a "mass term" to the scalar sector of gravity and modifies gravi- 
tational force in the linearized level even in Minkowski and de Sitter spacetimes. The 
Higgs phase of gravity provided by the ghost condensation is simplest in the sense 
that the number of Nambu-Goldstone bosons associated with spontaneous Lorentz 
breaking is just one and that only the scalar sector is essentially modified. 

In ghost condensation the linearized gravitational potential is modified at the 
length scale r c in the time scale t c , where r c and t c are related to the scale of sponta- 
neous Lorentz breaking M as 



Note that r c and t c are much longer than 1/M. The way gravity is modified is peculiar. 
At the time when a gravitational source is turned on, the potential is exactly the same 
as that in general relativity. After that, however, the standard form of the potential 
is modulated with oscillation in space and with exponential growth in time. This is 
an analogue of Jeans instability, but unlike the usual Jeans instability, it persists in 
the linearized level even in Minkowski background. The length scale r c and the time 
scale t c above are for the oscillation and the exponential growth, respectively. At the 
time ~ t c , the modification part of the linear potential will have an appreciable peak 
only at the distance ~ r c . At larger distances, it will take more time for excitations 
of the Nambu-Goldstone boson to propagate from the source and to modify the 
gravitational potential. At shorter distances, the modification is smaller than at 
the peak position because of the spatial oscillation with the boundary condition at 
the origin. The behavior explained here applies to Minkowski background, but in 

^ee e.g.[13|211|23EniEIllS2E3I for other related proposals. 

2 Lorentz violation in particle physics has been an active field of research. See e.g. 1551 1351 1771 
I38j and references therein. Gravity in the presence of a vector field with fixed norm is also extensively 
studied [3^1 EDI EH E21 E3] • For relation between this theory and the gauged ghost condensation, see 
Appendix of ref. |55). 
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ref. j2H] the modification of gravity in de Sitter spacetime was also analyzed. It was 
shown that the growing mode of the linear gravitational potential disappears when 
the Hubble expansion rate exceeds a critical value H c ~ l/t c . Thus, the onset of the 
IR modification starts at the time when the Hubble expansion rate becomes as low 
as H c . 

If we take the M/M-p\ — > limit then the Higgs sector is completely decoupled from 
the gravity and the matter sectors and, thus, the general relativity is safely recovered. 
Therefore, cosmological and astrophysical considerations in general do not set a lower 
bound on the scale M of spontaneous Lorentz breaking, but provide upper bounds on 
M. If we trusted the linear approximation for all gravitational sources for all times 
then the requirement H c < H would give the bound M < (M^Hq) 1 ^ 3 ~ IQMeV, 
where Hq is the Hubble parameter today [2E]. However, for virtually all interesting 
gravitational sources the nonlinear dynamics dominates in time scales shorter than the 
age of the universe. As a result the nonlinear dynamics cuts off the Jeans instability 
of the linear theory, and allows M < lOOGeV [44J. 

Many other aspects of ghost condensation have been explored. They include 
a new spin-dependent force [35], a qualitatively different picture of inflationary de 
Sitter phase jlHlllZj, effects of moving sources jlHlESl) nonlinear dynamics |oTH Hlj. 
properties of black holes [!)T1E21EB], implications to galaxy rotation curves |531I55| I5S]. 
dark energy models [S3 EH EH] , other classical dynamics jnOJEI], attempts towards 
UV completion [B21 EH] , and so on. 

In the simplest setup of the ghost condensation an exact shift symmetry is as- 
sumed and there is no potential term in the Lagrangian of the Higgs sector. As a 
result the Higgs sector behaves like a cosmological constant plus cold dark matter 
for homogeneous, isotropic background evolution 04]. If the shift symmetry is not 
exact but is softly broken then a shallow potential is allowed in the Higgs sector La- 
grangian. Recently, Creminelli, et.al [64J showed that the ghost condensation with 
softly broken shift symmetry can violate the null energy condition without any insta- 
bilities 3 . This opens up interesting possibilities of non-standard cosmology, including 
an accelerating universe with w < — 1. 

A point is that the coefficient of the time kinetic term 7r 2 in the low energy effective 
action for the scalar excitation 7r is positive and of order unity. This means that there 
is no ghost in the ghost condensation. This also implies that in ghost condensation 
there is no problem analogous to the strong coupling issues which the massive gravity 
theory and the DGP brane model are facing with. On the other hand, the coefficient 

3 Another model which violates the null energy condition without UV instabilities is proposed in 
|65| . It is interesting to notice that this model also breaks Lorentz symmetry. 
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of the space kinetic term (Vtt) 2 in the action becomes positive when w < — 1. (A 
usual non-ghost scalar field, with w > — 1, has a positive coefficient for the time 
kinetic term and a negative coefficient for the space kinetic term.) This does not 
necessarily introduce instabilities since there is also a higher-derivative space kinetic 
term (V 2 7r) 2 /M 2 with a negative coefficient. Actually, as far as the positive coefficient 
of (Vtt) 2 is small enough, the higher-derivative space kinetic term pushes the would- 
be unstable modes outside the cosmological horizon so that any instabilities do not 
show up. This situation is realized if the violation of the null energy condition is 
not too large since the positive coefficient of (Vtt) 2 is proportional to the amount of 
violation of the null energy condition [64J. 

In the present paper we investigate the classical dynamics of cosmology in ghost 
condensation with softly broken shift symmetry in more detail. Throughout this 
paper we shall adopt a AD covariant action explained in Sec. El In Sec. 01 we show 
that it is always possible to find a form of the potential in the Higgs sector Lagrangian 
which realizes an arbitrary FRW cosmological history {H(t), p(t)). A similar result 
is known to hold also for a conventional scalar field with a potential, but in this case 
it is impossible to violate the null energy condition without a ghost. This is the 
origin of the folklore that the phantom (w < —1) cosmology requires a ghost, which 
is correct for ordinary scalars. On the other hand, for the ghost condensate, the null 
energy condition can be violated without introducing ghosts or any other instabilities 
as far as the violation is weak enough [64J. Thus, the folklore is not correct in ghost 
condensation. It is probably worth stressing here again that the low energy effective 
field theory of ghost condensation is completely determined by the symmetry breaking 
pattern and does not include any ghosts. 

After showing the reconstruction method in Sec. El we investigate cosmological 
perturbation around the FRW background in Sec.El The resulting evolution equation 
is summarized in subsection 14.41 and can be used to test the theory by observational 
data of e.g. cosmic microwave background anisotropy, weak gravitational lensing and 
galaxy clustering and so on. 

2 AD covariant action 

Since the structure of the low energy effective field theory of ghost condensation is 
determined by the symmetry breaking pattern, it is not compulsory to consider a AD 
covariant description for the ghost condensation as far as physics in the Higgs phase 
is concerned. Nonetheless, it is instructive and sometimes convenient to have such a 
description. In this paper we shall adopt the AD covariant description and extend the 
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low energy field theory developed in ref. [22] to a general FRW background driven 
by not only the ghost condensate itself but also other cosmo logical fluids 4 . In this 
language, we start with a 4D covariant action principle for a scalar field, and the 
ghost condensation is realized as a background with a non-vanishing derivative of the 
scalar field which does not vanish even in Minkowski or de Sitter spacetime. Then the 
low energy effective field theory is obtained by expanding the covariant action around 
this background. Needless to say, the two approaches, one based on the symmetry 
breaking pattern and the other based on the 4D covariant action, completely agree 
when applied to physics within the regime of validity of the effective field theory. 

The leading 4D covariant action for ghost condensation is given by 
"M 4 



where 
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M is the symmetry breaking scale giving the cutoff scale of the low energy effective 
theory, and the sign convention for the metric is ( — h ++)• Note that the potential 
M A V((j)) is included to take into account the soft breaking of the shift symmetry |fi4j : 
V(4>) is assumed to depend on very weakly. When V(4>) is a constant, the shift 
symmetry is exact. The corresponding stress-energy tensor is easily calculated as 
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4 See Sec. 3.3 of Ref. [HI] for the low energy effective field theory in a FRW background driven 
solely by the ghost condensate itself. In the present paper we include general gravitational sources as 
well as the ghost condensate since the inclusion of those sources is essential for the test of the theory 
by e.g. cosmic microwave background anisotropy, weak gravitational lensing, galaxy clustering and 
so on. 
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where is the stress energy tensor of the other gravitational sources satisfying the 
conservation equation V V T" = 0. As a consistency check, it is easy to confirm that 

V V TW" = E^d^, (2.6) 

where is the left hand side of (|2.4jl . Thus, the stress-energy tensor of 4> satisfies 
the conservation equation V v Tjf^ v = 0, provided that the equation of motion i£w = o 
is satisfied. On the other hand, if the Einstein equation is satisfied and if d^cf) is non- 
vanishing then the equation of motion E^> = follows. Throughout this paper we 
shall use these expressions for the stress-energy tensor and the equation of motion. 

If we set «i = «2 = then the model is reduced to a kind of k-inflation [HE] 
or k-essence |6T1 IHHl ESI- However, in this case there is no modification of gravity 
in Minkowski or de Sitter spacetime. Moreover, with a% = ct2 = 0, the attractor 
£ = is unstable against inhomogeneous perturbations. In the presence of the 
terms proportional to a\ and a>2 {ol\ + ct2 > 0), the attractor is stable against small 
perturbations and gravity is modified in the linearized level even in Minkowski and de 
Sitter backgrounds. The relation between the k-inflation and the ghost condensation 
is in some sense similar to that between the usual potentially-dominated inflation and 
the Higgs mechanism [70] . 



3 FRW background 



Ghost condensation provides the simplest Higgs phase of gravity in which there is only 
one Nambu-Goldstone boson associated with spontaneous Lorentz breaking. Note 
again that the dynamics of the Higgs phase of gravity in ghost condensation has 
nothing to do with ghosts. Indeed, there is no ghost within the regime of validity 
of the effective field theory. Moreover, the Higgs phase of gravity has universal low 
energy description determined solely by the symmetry breaking pattern. 

In this section we consider a flat FRW ansatz 

g^dx^dx 11 = —dt 2 + a{t) 2 8ijdx l dx\ 

= <P(t), (3.1) 

where i — 1,2, 3, and analyze the dynamics of the homogeneous, isotropic universe. 
Linear perturbation around this background will be considered in the next section. 
With this ansatz, the stress energy tensor for the field <fi is 
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Here, we have expressed d t (f>, d 2 (p and (9 t 3 in terms of E = (<9 t 0) 2 /M 4 — 1 and its 
derivatives. The Einstein equation is 



3M 2 y H 2 
2M^d t H 



P<l> + P, 
-(pa + 



- (P+P), 



(3.5) 



where H = dta/a, and p and p are the energy density and the pressure of the other 
gravitational sources. Following the usual convention we call the first equation the 
Friedmann equation and the second the dynamical equation. Note that, because of 
the identity (j2.6|) and the Bianchi identity, the equation of motion for automatically 
follows from ()3.5|) . provided that the conservation equation for p and p holds: 



d t p + 3H(p + p) = 0. 



(3.6) 



3.1 Low energy expansion 

At this point, one might think that the expressions in (j3.3J) are too complicated to 
extract the physical picture of cosmology in the Higgs phase of gravity. Even if one 
could somehow manage to do so with brute force for the FRW background, it might 
be too optimistic to expect that the same approach works for linear perturbation 
around the FRW background. So let us go back and reconsider the meaning of the 
covariant action (j2.1|) . This is not a full action including a UV completion but just 
a leading action suitable to describe physics sufficiently below the cutoff scale M 
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around backgrounds with S ~ 0. If the Hubble expansion rate and/or the stress- 
energy tensor of the field <fi become close to or above unity in the unit of M then the 
low energy description is invalidated and we need a UV completion. Needless to say, 
the same criterion applies to the other approach based on the symmetry breaking 
pattern, and the two approaches agree in the regime of validity of the low energy 
effective theory. Therefore, all we can and should trust is what is obtained below the 
cutoff scale. In other words, we assume the existence of a good UV completion but 
never use its properties. For this reason, we can and should ignore terms irrelevant 
at low energies compared with the cutoff scale M. 

To be systematic, we adopt a low energy expansion by introducing small dimen- 
sionless parameters (i = 0, 1, 2, 3) as 

M 2 H p 2 p _ 2 

M|; = e °' M =£l ' Mp^ = £2 ' Mp^ = 0(e2) ' (3 ' 7) 

and 

Max(|S|,|V|) = e 3 , (E = 0(e 3 ), V = 0(e 3 )) . (3.8) 

Since the time scale for the change of H, E, V, p and p is expected to be the cosmo- 
logical time scale 1/H, it also follows that 

r) n H f) n T f) n V 

^l = 0^ + \ ^ = OW), ^ = OW)> (3-9) 



and that 



for n = 1, 2, • • •. The latter condition (J3.1U|) is not necessary for the present purpose 
but will be used when we derive evolution equations for linear perturbation around 
the FRW background in the next section. Unless fine-tuned, consistency of this 
assignment with the Einstein equation ()3.5|) requires that 

e 2 ~Max(e e 3 ,e 2 ). (3.11) 

It is in principle possible but not practical to perform the low energy expansion 
with respect to all e's, considering each as independent small parameters. Of 
course, if one performs the low energy expansion up to sufficiently high order with 
respect to all e's then one can always reach the point where all relevant terms (and 
probably many other irrelevant terms) are included. However, this is not the most 
economical way to obtain results relevant to particlar situations of physical interest. 
It is more economical and convenient to suppose some rough relations among these 
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small parameters to reduce the number of independent small parameters. This also 
makes it easier to extract physical picture out of complicated equations. Note that 
those relations among ej must reflect the situations of physical interest. 

For example, in ghost inflation we set £ = V = and include another field, 
say if), to end the inflation and to reheat the universe a la hybrid inflation. In this 
case e 3 vanishes. Moreover, since the stress energy tensor is dominated by the other 
field if), the consistency relation (|3.11|) implies that e± ~ e 2 . In this way, in Ref. [IE] 
we considered eo and e\ (~ e 2 ) as two independent small parameters. 

On the other hand, in the present paper we would like to consider late time 
cosmology in which both (p^,, p^) and (p, p) may contribute to the background as 
gravitational sources. For this reason and because of the consistency relation Ij3.ll}) . 
we suppose that there is a small number e such that 

4 = 0(e 2 ), e e 3 = 0(e 2 ), ej = 0(e 2 ), (e.g. e = Max( v ^, e 2 )) . (3.12) 

We expect (and will actually show) that in the dispersion relation for the excitation 
of ghost condensation, M 2 /Mp L (= e ) and S (= e 3 ) additively contribute to the 
coefficient of the momentum squared. (See the expression for Cq in 1)4.16}) below.) 
We would like to consider situations in which they can independently become rele- 
vant since they may become different sources of possible instabilities and interesting 
physical effects jHlj. Thus, we refine the second relation in ()3.12}) to 

e = O(e), e 3 = 0(e), (e.g. e = Max(e , e 2 , e 3 )) . (3.13) 

In summary we can introduce just one small parameter e and suppose that 

^L = 0(e), AiL = 0(e" +1 ), %fl = 0(6" +1 ), ^- = 0(e n+1 ), (3.14) 
M|j M n+1 M n M n 



and 



-0(e n+1 ), J^ = 0(0, (3-15) 



M n + A v ' M n+A 
where n = 0, 1, 2, • • •. Again, the latter condition ([3.15}) is not necessary for the present 
purpose but will be used when we derive evolution equations for linear perturbation 
around the FRW background in the next section. 

The low energy expansion not only enables us to compute various quantities in a 
systematic way but also helps us avoid picking up spurious modes associated with the 
higher derivative terms. In the 4D covariant action 1)2.1}) the terms proportional to 
a% and a 2 include the square of the second time derivative of the field 0. Therefore, 
if we take its face value then the equation of motion for includes up to the fourth 
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order time derivatives and there is in principle freedom to specify <fr, dt4>, d^cf) and d^cf) 
as an initial condition. However, as the scaling analysis in ref. shows, the time 
derivatives higher than the second order in the equation of motion are irrelevant at 
energies sufficiently below the cutoff M, at least in Minkowski background. In other 
words, extra modes associated with those higher time derivatives have frequencies 
above the cutoff scale M and is outside the regime of validity of the effective field 
theory. For this reason, those extra modes are spurious and must be dropped out from 
physical spectrum of the low energy theory. In the expanding background it is not a 
priori completely clear how to drop the spurious modes while maintaining all physical 
modes. By adopting the low energy expansion this can be done in a systematic way. 

3.2 Reconstructing the potential from H(t) 

As shown in Ref. [fH], many non-standard cosmology, including the phantom (w < 
— 1) cosmology, can be realized in the framework of ghost condensation without in- 
troducing ghosts or any other instabilities. The purpose of this subsection is to show 
that, given an arbitrary history of the Hubble expansion rate H(t) and gravitational 
sources pit) and pit), it is indeed possible to find a form of the potential V(<f>) which 
realizes H(t) as a solution to the Einstein equation (|3.5|) . provided that the conserva- 
tion equation d t p + 3H(p + p) = is satisfied. As is well known, a similar result holds 
also for a conventional scalar field with a potential: one can almost always find a 
form of the potential for a given history of the Hubble expansion rate. However, for a 
conventional scalar field, it is impossible to violate the null energy condition without 
a ghost. This is the origin of the folklore that the phantom (w < —1) cosmology re- 
quires a ghost, which is correct for ordinary scalars. On the other hand, for the ghost 
condensate, the null energy condition can be violated without introducing ghosts or 
any other instabilities as far as the violation is weak enough. Thus, the folklore is 
not correct in ghost condensation. It is probably worth stressing here again that the 
low energy effective field theory of ghost condensation is completely determined by 
the symmetry breaking pattern and does not include any ghosts. 

For the reconstruction of the potential V(4>) from the expansion history H(t) and 
the energy density p(t) and pressure p(t) of known gravitational sources, we take 
advantage of the low energy expansion introduced in the previous subsection. We 
expand V and E with respect to e: 

V = V + V! + ---, V n = 0(e n+1 ), 

E = E + E 1 + --- ) E n = 0(e n+1 ). (3.16) 
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Note that <j> is expressed in terms of S n as 

= M 2 J Vl + Edt = M 2 J 



1 + -S + 0(e 2 



dt. 



(3.17) 



The Friedmann and dynamical equations ()3.5|) in the lowest order in e are easily 
solved with respect to S and Vq as 



So 
V 



M 4 



2M 2 1 <9 t #+(p + p) 



M 4 



M^(2d t H + 3H 2 ) + p 



The relation (|3.17jl in the lowest order in e is 

(f) = M 2 (t-t ) + M-0(e), 



(3.18) 



(3.19) 



where t is a constant. Thus, for a given history (H(t), p(t), pit)) satisfying d t p + 
3H(p +p) — 0, the potential V is expressed in terms of (f) as 



V(<j>) = ^ I [M 2 l (2d t H + W 2 )+p 



t=t +4>/M< 



+ 0(e 2 



(3.20) 



It is also easy to obtain the correction to the lowest order potential. In the next- 
to-the-leading order the Friedmann and dynamical equations f!3.5|) are solved with 
respect to Si and V\ as 



Si 



-^-w dtH+ J^ {5dtH+3H2) 



1 o 3a — 23 , > 



2M 2 



-(2d t H + 3H 2 ). 



The relation (|3.17|) in this order is 



M 2 rt 

M 2 {t - t ) + — / X (t')dt' + M ■ 0{e 2 
2 Jt a 



(3.21) 



(3.22) 



where to is again a constant. Thus, for a given history (H(t), p(t), p{t)) satisfying 
d t p + 3H(p + p) — 0, the potential V is expressed in terms of <fi as 



V((f>) = [V + At ■ d t V + VI] 



+ 0(e 3 



where 



At 



2 A, 



E (t')dt' 



(3.23) 
(3.24) 



and S , Vq and Vi are given in (|3.18|) and (|3.21j) . 
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4 Cosmological perturbation 



In this section we derive evolution equations for cosmological perturbation around 
the FRW background. 

We consider a scalar-type cosmological perturbation in the longitudinal gauge, 



ds 2 = -(1 + 2^Y)dt 2 + a 2 (l + 2$y)dx 2 , 
and a general stress-energy tensor of the form 



(4.1) 
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(4.2) 



where F's are harmonics on the 3-space defined as 
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Yi 
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ay, Y i = 6 ij Y, 
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5 ik Y kj . 



(4.3) 



Here, p and p are the unperturbed energy density and pressure, and (A, V, T, II) rep- 
resent the gauge-invariant perturbation of the stress-energy tensor. Physical meaning 
of each component is as follows [21]: A is the density contrast in the slicing that 
the fluid velocity is orthogonal to constant time hypersurfaces, V is the fluid veloc- 
ity relative to the observers normal to constant time hypersurfaces, T is the entropy 
perturbation and II is anisotropic stress. They satisfy the perturbed conservation 
equation 



d t (pA) + 3HpA + — - 3d t H 



(p + p)V + 2H P U 



+3(p + p)(dt$ - HV) 



dt [(p + p)V] + (4 + 3c 2 )H(p + p)V 



pT + c 2 s pA + (p + p)W - -pU 



0. 

0. 

(4.4) 



When p and p + p are non- vanishing, these equations are rewritten as 

k 2 



d t A-3wHA + (1 + w) 



3d t H 



V + 3{d t $ - Hty) 



+ 2wHU 



12 



dtV + HV-^ 



1 + w 1 + w V 3 



0. 

(4.5) 



where w = p/p and cj. = d t p/d t p. 

As for the metric perturbation, for convenience, we decompose into the stan- 

dard, general relativity (GR) part (&gr^gr) an d the modification part {$ mo d^mod) 
as follows. 

$ = $ G R + $mod, 

V = ^GR+^mod, (4.6) 



where the GR part is given by 



_pA_ 

a 2 VGR , 



'PI 

,2 



k /T ^ N P U 



-(*gr + <S>cr) = - W . (4.7) 



In the following we shall derive evolution equation of the modification part of the 
metric perturbation (& mo d, ^mod)- For readers who are interested in application of 
the formalism, subsection 14.41 summarizes resulting evolution equations of the cos- 
mological perturbation in the Higgs phase of gravity. In subsection 14.11 we introduce 
a systematic low energy expansion applicable to linear perturbation by extending 
the low energy expansion for the background given in subsection 13.11 Subsection 14.21 
includes the hardest part of calculations in this paper: we derive a single master 
equation governing the modification of gravity without ignoring higher time deriva- 
tives. The master equation is a fourth-order ordinary differential equation for each 
comoving momentum. However, in subsection 14.31 we see that the fourth and third 
order derivative terms are irrelevant and reduce the master equation to a set of two 
first-order ordinary differential equations for each comoving momentum. At the same 
time, we remove an apparent singularity in the master equation. In subsection 14.41 
we summarize the results of this section in a way which is directly applicable to ac- 
tual problems. In Appendix IA.1[ as a simple application of the formula summarized 
in subsection 14.41 we consider Minkowski and de Sitter backgrounds and derive the 
results in ref. for modification of gravity in those backgrounds. Readers who are 
interested in application of the formalism may go directly to subsections 14.41 
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4-1 Low energy expansion 



We now extend the low energy expansion developed in subsection 13.11 to linear per- 
turbation around the FRW background. The assignment for the FRW background is 
summarized in (|3.14j) and (|3.15j) . 

For perturbation, besides the smallness of the perturbation itself controlling the 
validity of the linear approximation, there is one additional small parameter e 4 given 
by 

1 , x 

M— = < 4 ' 8 > 
The low energy effective theory is valid only if e 4 is small enough. Supposing that 
A, V, T and IT may become comparable, the perturbed conservation equation (j4.5|) 
implies that 

^(A, v, r, n) ~ (A, v, r, n) • Max( ei , e 4 ), (4.9) 

unless the modification part ($ mo d, ^mod) dominates over the GR part (&gr, *&gr)- 
As shown already in ref. j2E], the dispersion relation for excitation of ghost conden- 
sation in Minkowski background is uo 2 ~ «k 4 /M 2 — ak 2 M 2 /2M| 1 . The first term in 
the right hand side is ~ e\M 2 and the second term is ~ ee 2 M 2 . In this paper we 
would like to generalize this dispersion relation in Minkowski spacetime to an evo- 
lution equation in the FRW background, taking both terms into account. For this 
purpose we suppose that these two terms may become comparable. In other words, 
we suppose that there is a small parameter e such that 

e = 0(e 2 ), e 4 = 0(e), (e.g. e = Max( v ^, e 4 )) . (4.10) 

In summary it is supposed that there is a small parameter e such that 
M 2 



Oil 



~2\ 



M 2 Pl ^ " 



r) n H f) n y f) n V 

U t n = n( ~2(n+l)\ <±±_ _ n( ~2(n+l)\ °t v _ n( ~2(n+iy 

M n+l - W ^ )> M n W ^ )> M n ^ 

®tP _ r^(~2{n+l)\ d?P _ n f~2(n+l) 



0(^ n+1) ), T^ = 0(e 2( " +1) ), (4-11) 



and 



M n + 4 ' M™+ 4 



0(1 




M a 

(A,v,r,n) ~ o(e")-(A,y,r,n), (4.12) 



where n — 0, 1, 2, 
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4-2 Master equation 



In the previous subsection we developed a low energy expansion for the FRW back- 
ground and linear perturbation around it. Armed with this, we can roughly estimate 
how small or large a term in equations should be, as far as (&gr, ^gr), (A, V, T, 
IT) and their time derivatives are concerned. On the other hand, we do not know a 
priori how small or large the time derivatives of the modification part (§ m od, ^mod) 
are. This can be seen only after we derive evolution equation of these quantities. In 
this section we derive a single fourth order master equation for the modification part, 
using the low energy expansion but not assuming the smallness of the time derivatives 

Of (®mod, ^mod)- 

Let us introduce two variables $1,2 by 



where c is a constant. We shall see below that the source terms for equations governing 
$1,2 scale like oc M 2 /M^ in the M/M P1 ->• limit. Thus the standard GR result is 
recovered in the M/Mp\ — > limit. The variable $2 (oc Q mo d + ^mod) was defined so 
that it vanishes when (3 = 0. The definition of the variable $1 is not yet fixed because 
of the unfixed constant c. The value of the constant c will later be determined so 
that the master equation for <3>i including matter source terms is simplified. 

As for the field <fi responsible for ghost condensation, we expand it up to the linear 
order as 



We now have three unknown variables $1, $2 and it other than matter variables (A, 
V, T, II). Our task now is to obtain the evolution equation for those three variables 
sourced by these matter variables. 

In cosmology with an ordinary scalar field the standard strategy to analyze lin- 
ear perturbation in the longitudinal gauge is as follows: (i) to eliminate one of the 
two metric variables ($, say ^, by using the traceless part of (ij)-components 
of the linearized Einstein equation; (ii) to eliminate the linear perturbation of the 
scalar field by using the (Oi)-components of the linearized Einstein equation; and (iii) 
to obtain a single master equation for the remaining metric variable, say $, from 
the (OO)-component of the linearized Einstein equation. (The linearized equation of 
motion of the scalar field and the trace part of (ij)-components of the linearized 
Einstein equation are automatically satisfied because of the Bianchi identity.) The 
steps (i) and (ii) involve simple algebraic equations for variables being eliminated. 





$1 + C$; 



2, 



(4.13) 




(4.14) 
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As a result the master equation is a second-order ordinary differential equation for 
each comoving momentum. This is consistent with the fact that, in the absence of 
other gravitational sources, the scalar sector includes only one propagating degree of 
freedom, i.e. excitation of the scalar field. If there are other gravitational sources 
then the master equation includes source terms due to those gravitational sources. 

We can use the same strategy to analyze the linear perturbation with ghost con- 
densation although the result will be drastically different. There are again three main 
steps. (To be precise, for the reason explained below, steps (ii) and (iii) are not sep- 
arated but actually mixed in the present case.) (i) First we eliminate $2 by using 
the traceless part of (^-components of the linearized Einstein equation, (ii) Next 
we eliminate it by using the (Oi)-components of the linearized Einstein equation, (iii) 
Finally we obtain a single master equation for the remaining metric variable $1 from 
the (OO)-component of the linearized Einstein equation. (The equation of motion of 
7r and the trace part of (^-components of the linearized Einstein equation are au- 
tomatically satisfied because of the Bianchi identity.) The step (i) involves just an 
algebraic equation for $2 as in the standard case. However, the step (ii) involves a 
second-order differential equation for n for each comoving momentum, contrary to 
the standard case. For this reason the steps (ii) and (iii) are not separated but ac- 
tually mixed. After eliminating i\ and its derivatives, the resulting master equation 
is a fourth-order ordinary differential equation for $1 for each comoving momentum. 
This is indeed consistent with the fact that the action for includes the square of the 
second derivative of (f> and that the equation of motion for n involves up to forth-order 
derivatives at least formally. In the master equation the matter variables (A, V, T, 
IT) appear as source terms. 

We adopt the low energy expansion introduced in the previous subsection to follow 
each step (i)-(iii). This in particular means that terms with sufficiently higher time 
derivatives acted on the background or matter variables can be ignored since they are 
considered as higher order in the expansion with respect to e. On the other hand, 
the time derivative acted on $1, $ 2 and n is not supposed to raise the order of the 
low energy expansion. This is because we do not a priori know the time scale of the 
dynamics of these three variables in the FRW background until the evolution equation 
of a single master variable is obtained. 

Following the steps (i)-(iii), we obtain the fourth-oder master equation for $!: 

+ C 3 d t 3 $! + + Ci^di + C $i = S u (4.15) 
where the coefficient of the second-order term is normalized to unity, and other coef- 
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_a_ Oje 2 ) 

M 2 + ^f 2 ~' 



ficients and the source term are given by 
C 4 
C 3 

Ci 
C 



a 

M 2 



6H- 



M 2 (d t Y> + 2HY 1 ) 
M 2 £ + 2ak 2 /a 2 



+ 



0(i 



M ' 



= 3H 



M 2 (<9 t £ + 2#£) 
M 2 £ + 2ak 2 /a 2 



+ M-0(e 4 ), 



a k 4 1 / M 2 \ k 2 rr2 a TT 



M 2 a 4 ' 2 



Pi/ 



SM 4 M 2 H(d t T, + 2HT) 



+ M 2 -0(e € ), 



4M 2 ! M 2 E + 2ak 2 /a 2 



and 



k 2 ^ 



PI „2 

or 



5i = 



M 2 

Mfi I 8 



i (M 2 S + 2ak 2 /a 2 ) [(1 + 3c 2 )pA + 3pr 



+ cPd 2 ( P U) 



+(c-l)/3 
+M 6 -0(e 5 ). 

On the other hand, $2 is expressed in terms of $1 and its derivatives as 
$ 2 = (D 3 d^! + D 2 d 2 $ 1 + D ± dA + A)^ + S 2 ) , 



where 



D 3 = 

D 2 = 

D 1 = 

D n = 



AaH 



M 2 [M 2 S + 2ak 2 /a 2 ] ' M 3 



+ 



2a 



M 2 [a - (2c + 1)13} ' M 2 ' 



+ 



AH 



+ 



0(e 2 ) 



M 2 E + 2ak 2 /a 2 M ' 
2 k 2 M 2 AH 2 

+ T7T- + 



M 2 l —S 2 



M 2 a 2 M 2 S + 2ak 2 /a 2 

M 2 



+ 0(e 



~4s 



2M 2 ! 



(pA + 2pn) + M 4 -0(r). 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



To obtain these expressions, we kept terms up to 0(e 4 ) in intermediate steps. 
Having obtained the master equation, it is clear that 



c= 1 



(4.20) 



17 



gives the best choice for the definition of $1. This indeed simplifies the source term Si 
of the master equation. Moreover, the term (M 2 / Mp t ) [3d 2 (pit) in Si can be absorbed 
by redefinition of the master variable: it does not appear if we use $1 and $2 defined 
by 

1 " 1+P k 2 M| 1 ' 
$2 = $2, (4.21) 

instead of $1 and $2- The master equation is now written as 

C 4 <9 t 4 $i + C 3 af$i + d^i + dd t ^i + C $i = Si, (4.22) 

where C^x,a are given by ()4.1fij) and the new source term is 

M^%Si = -^r (M 2 E + 2«k 2 /a 2 ) [(1 + 3c 2 )pA + 3pT] + M 6 • 0(e 5 ). (4.23) 

The variable $2 is expressed in terms of $1 and matter variables as 

$2 = P (D^i + D 2 d^i + Di^i + A)$i + 5a) , (4.24) 

where D 3 2 ,i,o an d 5" 2 are given by (|4.19|) with c = 1. 

Note that coefficients of the master equation ()4.22j) and the relation (|4.24j) were 
obtained just up to the leading order in the e expansion. It is of course possible to seek 
corrections to the leading terms, but we shall not do so in this paper for simplicity. 
Note also that we had to keep up to the fourth order in e in the intermediate steps 
to obtain the leading master equation. The reason for this is manifest: otherwise, we 
would not have been able to obtain a non-vanishing coefficient of $1 in ()4.22|) since 
Co = 0(e 4 ) while the coefficient of <9 2 $i is unity. 



4-3 Reduction to a set of two first order equations 

Although the master equation ()4.22|) obtained in the previous subsection is formally 
fourth order, the fourth and third order derivative terms are actually irrelevant. To see 
this is easy: they are suppressed by the cutoff scale M compared with the coefficients 
of the second and first derivative terms, respectively. This means that ignoring the 
fourth and third order derivative terms just drop out modes with frequency of order 
M or higher, which are outside the regime of validity of the effective field theory. 
Properties of these modes are dependent of properties of unknown UV completion 
and, thus, we not only may but also must discard these modes. In other words, 
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we assume the existence of a good UV completion but will never use its properties. 
Hence, we obtain the second-order master equation 



(4.25) 



where Ci t o and Si are given by ()4.1fi|) and ()4.23j) . For the same reason, we must 
discard the third and second order derivative terms in the relation ()4.24|) : 



(4.26) 



where -Di j0 and S 2 are given by (j4.19j) with c = 1. 

The master equation (J4.25)) and the relation (J4.26|) appear to be singular when 
M 2 £ + 2ak 2 /a 2 vanishes. This does not mean that the system exits the regime 
of validity of the low energy effective theory but just implies that the dynamics 
is not described by a single master equation. Indeed, at the would-be singularity 
M 2 £ + 2ak 2 /a 2 = 0, the second-order master equation is reduced to d t $i +H$i = 0. 
This suggest that the second-order master equation should be decomposed into a set 
of two first-order equations. To make this more explicit, let us define a new variable 

X M 2 E + 2ak 2 /a 2 ' 1 ' ; 

With this variable and $i the master equation ()4.25|) is reduced to the following set 
of two first-order equations. 



dtX = 

where the source term S x is given by 



M 2 S + 2a- 



M 2 

Mi 



_2_kj^ 

M 2 a 2 ) 



$i + S- 



M^S X 



M 2 
2M 2 ! 



1 + 3c 2 )pA + 3pT 



On the other hand, the variable $2 is written in terms of $1 and \ as 

M 2 2 k 2 ^ 



M 2 ! 



M 2 



$1 + S 2 



where 



k 2 

M 2 ,-S 2 



M 2 
2M 2 l 



(pA + 2 P n). 



(4.28) 



(4.29) 



(4.30) 



(4.31) 
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The apparent singularity at M 2 £ + 2ak 2 /a 2 = does not exist in the new equations 
(|4.28|) and (|4.30|) . As easily seen from the structure of the set of first order equations 
(|4.28|) . the apparent singularity in (|4.25|) and (|4.26|) was just the reflection of the 
fact that the evolution of the variable $1 is decoupled from that of \ when M 2 £ + 



As shown in subsection 13 .2^ in ghost condensation it is always possible to find a form 
of a potential for the Higgs sector which realizes a given cosmological history (H(t), 
p(t), p(t)) of the FRW background, provided that the conservation d t p+3H(p+p) = 
is satisfied. An important point is that the null energy condition can be violated in 
the Higgs sector of gravity without introducing ghosts or any other instabilities as 
far as the violation is weak enough [U^. Note again that the structure of the low 
energy effective field theory of ghost condensation is completely determined by the 
symmetry breaking pattern and does not include any ghosts. Thus, for example we 
can realize the phantom (w < —1) cosmology without a ghost. This is an explicit 
counter-example against the folklore that the phantom (w < —1) cosmology would 
require a ghost. 

In this subsection we summarize the result of this section, where we have derived 
a set of two first order ordinary equations governing cosmological perturbations of the 
Higgs sector of gravity. The low energy effective theory of ghost condensation includes 
two dimensionless parameters of order unity a (= a± + «2 > 0) and (3 (= o^) and the 
cutoff scale M. Provided that (a, (3, M) are fixed and a background cosmological 
history (H(t), pit), pit)) (satisfying the conservation equation but possibly including 
deviation from general relativity) is given, the evolution equations for cosmological 
perturbations are completely determined. 

Now let us start summarizing evolution equations for cosmological perturbation. 
For notation the reader should refer to the second and the third paragraphs of this 
section. 

The metric perturbation is decomposed into the standard, general relativity (GR) 
part (§gr,^gr) an d the modification part {$ mo d^mod) as f|4.6|) . and the modification 
part is written as 



2ak 2 /a 2 



0. 



4-4 Results summarized 



mod 



$i + $ 2 - (3 



a 2 pll 




^ mod 



-2$ 2 +/3 



(4.32) 
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The variable $1 is given by solving 
dt^ + HQx = 

9tX = (^-^J*i + S x , (4.33) 
where \ is an auxiliary variable, S is given by 

and represents deviation of the background evolution from general relativity, and the 
source term S x is given by 




2M* 1 d t H +(p + p) (4.34) 



M pK S x = ^t[( 1 + PA + 3 P r] . (4.35) 



l a 2 * 2M P1 

On the other hand, the variable $2 is written in terms of $1 and x as 



$ 2 =/3 

where 



rr / M 2 2 k 2 \ „' 



(4.36) 



k 2 M 2 

M ^ = ^ A + W (4 ' 37) 
Note that $2 vanishes if (3 = 0. 

The source terms S x and S2 vanish in the M/Mp\ — > limit. This means that 
the modification part ($ mo d, ^mod) is not induced by matter sources in this limit. In 
other words, in this limit the Higgs sector of gravity is decoupled from the gravity 
and the matter sectors, and the general relativity is safely recovered. 

One must be aware that these equations are valid only if the physical momentum 
\f\s?/a is sufficiently lower than the scale of the spontaneous Lorentz breaking M, 
which plays the role of the cutoff scale of the low energy effective theory. 

The coefficient So includes information about deviation of the background evo- 
lution from general relativity. In usual approach, the deviation is considered as in- 
dication of dark energy and/or dark matter. Instead, we here replace those dark 
components with the Higgs mechanism of gravity, i.e. ghost condensation. The cases 
with S > and S < correspond to dark components with w > — 1 and w < — 1, 
respectively. In the usual approach, a dark component with w < — 1 is called phan- 
tom and is thought to be associated with ghosts, i.e. excitations with wrong-sign 
kinetic energy. On the other hand, as explained in Sec. ^ there is no ghost in the 
ghost condensation. The set of equations summarized above can be applied to both 
S > (w > —1) and S < (w < —1) cases, including transition between these two 
regimes. 
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5 Concluding remark 



In the simplest setup of the ghost condensation the shift symmetry is exact and it 
behaves like a cosmological constant plus cold dark matter for homogeneous, isotropic 
background evolution. With soft breaking of the shift symmetry, a shallow potential 
is allowed in the Higgs sector Lagrangian. We have investigated the classical dynamics 
of cosmology in ghost condensation with softly broken shift symmetry. As we have 
shown explicitly, it is always possible to find a form of the potential in the Higgs 
sector Lagrangian which realizes an arbitrary cosmological history (H(t), p(t)) of the 
visible sector of the FRW background. After showing the reconstruction method, we 
have derived the evolution equation for cosmological perturbations in the Higgs phase 
of gravity. 

The strongest evidence for the accelerating expansion of the universe today comes 
from the supernova distance-redshift relation. Using this kind of geometrical informa- 
tion of the universe, it is in principle possible to reconstruct the potential for the Higgs 
sector by the method developed in subsection 13.21 of this paper. Once this is done, 
the theory acquires predictive power. By using the formalism of the cosmological 
perturbations summarized in subsection 14 A\ this theory can be tested by dynam- 
ical information of large scale structure in the universe such as cosmic microwave 
background anisotropy, weak gravitational lensing and galaxy clustering. 

Results in this paper have been obtained within the regime of validity of the 
effective field theory, whose structure can be determined solely by the symmetry 
breaking pattern. Therefore, while we have used a covariant 4D action of a scalar 
field as a tool for calculation, the final results summarized in subsection 14.41 should 
be universal and independent of the way the Higgs phase of gravity is realized. The 
same equations should hold as far as the symmetry breaking pattern is the same, with 
or without a scalar field. 
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Appendix 



A.l Simple example: Minkowski and de Sitter backgrounds 

In this appendix, as a simple application of the formula summarized in subsection 14.41 
we consider modification of gravity in Minkowski and de Sitter backgrounds. This is 
the situation considered in ref. |26| . 

For simplicity we set (3 = 0. In this case $2 = and the modification part of the 
metric perturbation is 

®mod = -^mod = $1- (A.l) 



By setting 
and 



0, H = H (= const.) (A.2) 



pA = 5p, c 2 = T = IT = 0, (A.3) 
the set of first order equations ()4.33|) is reduced to 

O ^ TT ^ "If 2 

Ot^mod + H ® mod = 

2 or 

( M 2 2 k 2 \ 

dtx = [M-w^r^ 3 *' {aa) 

where the source term is now given by 

M 2 a 2 M 2 
S -=2MJ^ 5P= K^ GR - (A ' 5) 
By eliminating x from these equations, we obtain 

orro,, fcxk 4 CtM 2 k 2 n rr2 \ ^ «M 2 , . 

* *- + afflA *- + - 2M^ + = 2Jii£ (A - 6) 

This is eq. (8.27) of ref. 5 and explicitly shows that in the M/Mp\ — > limit, 
the GR part ceases to act as the source of the modification part $ mo d an d the 
general relativity is safely recovered. By introducing the length and time scales r c 
and t c as 

5 Note that a 2 in |2fi| corresponds to a in the present paper. 
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This equation is rewritten as 

/ 4 U4 2 1^2 \ 2 U2 

Vt 2 . a 4 t= a 2 J f 2 a 2 

In the Minkowski spacetime, by setting H = and a = 1 in (JA.8j) . it is easily 
seen that modes with the length scale ~ r c are unstable and that the time scale of 
the instability is ~ t c . This is the analogue of the Jeans instability found in ref. [2E] 
and explained in Sec. ^ What is different from the usual Jeans instability is that this 
behavior in the linearized level persists even in Minkowski spacetime 6 and modifies 
the linearized gravitational potential. 

It is also easy to see from (jA.8|) that the Jeans instability disappears when the 
Hubble expansion rate Ho is larger than a critical value H c ~ l/t c . Thus, the onset 
of the IR modification starts at the time when the Hubble expansion rate becomes 
as low as H c . 
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